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We present a fully analytically solvable family of models with many-body cluster interaction and Ising in¬ 
teraction. This family exhibits two phases, dubbed cluster and Ising phases, respectively. The critical point 
turns out to be independent of the cluster size n -|- 2 and is reached exactly when both interactions are equally 
weighted. For even n we prove that the cluster phase corresponds to a nematic ordered phase and in the case of 
odd n to a symmetry protected topological ordered phase. Though complex, we are able to quantify the multi¬ 
particle entanglement content of neighboring spins. We prove that there exists no bipartite or, in more detail, 
no n -I- 1-partite entanglement. This is possible since the non-trivial symmetries of the Hamiltonian restrict the 
state space. Indeed, only if the Ising interaction is strong enough (local) genuine n -|- 2-partite entanglement 
is built up. Due to their analytically solvableness the n-cluster-Ising models serve as a prototype for studying 
non trivial-spin orderings and due to their peculiar entanglement properties they serve as a potential reference 
system for the performance of quantum information tasks. 
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I. INTRODUCTION 


In many-body systems described by classical mechanics the 
presence of an ordered phase is connected to the spontaneous 
breaking of symmetries associated with local order parame¬ 
ters. A system consisting of classical spins, for instance, may 
admit a ground state having all spins oriented along a given 
direction. Such ground states simultaneously break the spin- 
rotation and the time-reversal symmetry witnessed by a non¬ 
vanishing magnetic moment. 

Considering quantum systems, in contrast, one finds also 
different phases connected to some physical quantity but not 
necessarily to the magnetic moment. The paradigmatic ex¬ 
ample is a translation invariant spin-i chain for which the 
ground states correspond to the so called valence bond states, 
i.e. states build up ^ tensor products of maximally entan¬ 
gled bipartite states Jlllj]- In such systems neither the spin- 
rotation nor the time-reversal symmetry is broken, neverthe¬ 
less, it is possible to define order parameters characterizing the 
phases. Typically nematic phases occur if at least one symme¬ 
try of the Hamiltonian is broken which are phases with long 
range ordering, i.e. defined by order parameters on a finite set 
of sides. Other exanmles intensively discussed are topologi¬ 
cal order phases IS 0] that, for instance, are associated with 
the robustness of ground state d^eneracies S] , are quantized 
non-Abelian geometric phases |0 or possess patterns due to 
long-range quantum entanglement 10 . 

Frustration occurs for systems with competing interactions 
or non-trivial geometries and can be related to quantum en¬ 
tanglement 101. Non-trivial spin orders are usually found if an 
interplay between frustration and quantum fluctuations is at 
work resulting in chiral, nematic or general multipolar phases. 
In contrast to topological phases, even in the case of a van¬ 
ishing magnetic moment the spin-rotation symmetry is bro¬ 
ken 10,0. These phases are also interesting from the point 
of applications. The topological ordered phases play a fun¬ 


damental role in the spi n liq uids |[iS[nl] and in non-Abelian 
fractional Hall systems lll^ and are predicted to play a key 
role in the future development of fault-tolerant quantum com¬ 
puters 10. The nematic order is usually found in materials 
commercially used in the liquid crystal technology 0 such 
as LCDs (liquid crystal display). 

Non-trivial ordered spins appear usually for higher dimen¬ 
sional systems (lattices) or sites with more than two degrees 
of freedom (spins higher than i). Exceptions are the frus¬ 
trated one dimensional ferromagnetic spin-i chain in an ex¬ 
ternal magnetic held having a nematic ordered phase 00 
and the one dimensional cluster-Ising model exhibiting a sym¬ 
metry protected topological ordered phase 00. 

In general, mathematical tools to handle such systems are 
rare and only few very specific Hamiltonians have been found 
to be analytically solvable. The present paper introduces 
a huge class of analytically solvable one dimensional mod¬ 
els with two degrees of freedoms (spin-i) exhibiting both 
topologically and nematic ordered phases, which we dub n- 
cluster-Ising models. The index n = 1,2,... refers to the 
presence of an n + 2-body interaction, a cluster size of n -F 2. 
The physical systems under investigation are characterized 
by two competing interactions, a two-body Ising interaction 
along the y-axis and an n H- 2-body interaction along the x- 
axis and the z-axis. The Hamiltonian of the family of models 
can be written as 


sin((^)^(TjcrJ+i -COs((^)^(TjO|_„crJ+„+i (1) 


where J has the dimension of an energy (which we set equal 
to one in the computation) and „ stands for 


n 





( 2 ) 
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Via the parameter cj) the relative weight of the two interactions 
is controlled; When (p approaches 0 the system is dominated 
by the multi-body interaction whereas when p tends to -k 12 
the system is dominated by the (anti-ferromagnetic) Ising in¬ 
teraction. 

We show that this family of models can be analytically 
solved (Sec.ini) and how the spin correlation function can be 
obtained (Sec. inii. We prove that there is a quantum critical 
point at (pc = 7r/4 separating the cluster phase from the Ising 
phase. This corresponds to the case when both interactions 
have equal weights. This critical point (pc, surprisingly, does 
not depend on the n -f 2-body interaction, hence, it does not 
dependon the cluster size. In strong contrast to the relevant or¬ 
dering in the cluster phase that depends strongly on n: In case 
of odd n a symmetry protected topologically ordered phase 
is present, whereas for even n a nematic phase is present. In 
both cases we determine the order parameter (string order pa¬ 
rameter for the topological ordered phase and block order pa¬ 
rameter for the nematic phase) as well as the order parameter 
of the Ising phase (Sec. 113 . 

In the next step we study the various entanglement proper¬ 
ties of the family of models (Sec. 0. The first observation is 
that for any n and (p - as proven for the standard cluster-Ising 
model (n = 1) in Ref. lU^ - there is no bipartite entangle¬ 
ment. Picking out any two spins the state is separable. In¬ 
deed, we find that this family of Hamiltonians lead to ground 
states that possess genuine k = n + 2-partite entanglement 
between any contiguous spins and any k < n + 2-partite en¬ 
tanglement vanishes. The symmetries in the state space of 
the ground states force the reduced state of any n + 2 adjacent 
spins into a so-called X-form ll^ . i.e. by applying certain lo¬ 
cal unitary operators the reduced density matrix has only non¬ 
zero entries on the two diagonals. Due to this form we can 
exactl y ev aluate a measure for genuine multipartite entangle¬ 
ment ll2ll - l^ . i.e. quantify the entanglement content. So far, 
long range multipartite entanglement close to a phase transi¬ 
tion has been studied in terms of entanglement witnesses, e.g . 
for the XXZ spin chain ll^ or for the XY model 12^2'^ . 
Having this strong tool at hand, a measure of genuine multi¬ 
partite entanglement, we find that non-zero genuine multipar¬ 
tite entanglement is only non-zero in the Ising phase (p > (pc 
(except n = 1), thus exhibiting a fortunate behaviour for ap¬ 
plications such as utilizing these quantum systems for quan¬ 
tum algorithms. 

The block entanglement properties are studied with focus 
around the quantum phase transition. Via the relation between 
conformal field theory ll^ and the divergence of the block 
entanglement at the quantum phase transition we are able to 
evaluate the central charges of the models that turns out to 
depend on n. 

Last but not least we conclude (Sec. IVIb by discussing the 
interplay between the characterization of the many-body sys¬ 
tems by ordered parameters and by symmetries in the Hilbert- 
Schmidt space of the ground states revealing the entanglement 
properties. 


II. SOLUTION OF THE MODELS 


In this section we present how to compute analytically the 
ground states of the models under investigation. The idea is 
to map the Hamiltonian, Eq. O, of spin-i particles into non¬ 
interacting fermions moving freely along the chain only obey¬ 
ing Pauli’s exclusion principle. This method works even for 
the case in which the length of the system diverges 
Having finally computed the energy density function we find 
a phase transition that is further analysed in Sec. |I3 

The mapping of a spin model to a fermionic one is obtained 
by applying the Jordan-Wigner transformation j^ . Provid¬ 
ing the correct anti-commutation rules in the Jordan-Wigner 
transformation one associates the local spin operators with 
non-local fermionic operators 


i-i i-i 



k=l k=l 


where are the respective ladder operators. Herewith the 
Hamiltonian in Eq. O becomes 

Jsin(0)^(cjc]+i - c]cj+i -f - c^Cj+i) 

3 

+ Jcos(0) CjCj+n+1 

3 

~ ^j'^j+n+1^ ( 4 ) 


One notes that herewith the cluster interaction is reduced from 
a n + 2 interaction to a two-body interaction between sites 
at distance n + 1. After having reduced the problem to an 
effective two-body one the model can be diagonalized via the 
Eourier transforms of the fermionic operators, i.e. 




-ikj 


Ck e 


d „ikj 


( 5 ) 


where the wave number k is equal to fc = 2ttI /N and I runs 
from —N/2 to N/2 and N is the total number of spins (sites) 
in the chain. The Hamiltonian transforms to 

= ^4") (6) 

fc >0 


with 

4 "^ = (blbl^-b.kbk) 

+2£fe,n ~ , 

where the parameters Sk,n and ek,n are respectively given by 

bk,n = J (sin ((n-b l)fc) cos(/)-f sin(A:) sin^) , 

£k,n = J (cos ((n-b l)fc) cos(^ — cos(fc) sin(()) . (7) 
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Via these transformations we re-wrote the Hamiltonian under 
investigation into the sum of non-interacting terms h!'^\ each 
one of them acting only on fermionic states with wave num¬ 
ber equal to k or —k. Each corresponds to a four level 
system that can be expressed in an occupation number basis 
by |lfe,l-fe), |0fc,0_fc), |lfe,0_fc), |0fc,l_fc) and is, explicitly, 
represented by the following matrix 

( £k,n 0 0\ 

7 (n) _ 2 i 6k^n 2 6k^n 0 0 I -Q. 

- 0 0 0 0 ’ 

Vo 000/ 



which ground state energy computes to 


l</> - 7r/4| 


^fc”^=-V^ln+'^L = “2'/\A-cos((n-f2)A:)sin(2(^). (9) 

The associated ground state ^ superposition of 

|lfc, l_fc) and |0fc,0_fc) 


Figure 1: (Color online) The graphes show the second derivative of 
the energy density of the ground state En, 4 , as a function of (f> for 
different cluster sizes n + 2. The divergence is independent of n at 
the critical value 4>c = ^ and corresponds to a vanishing energy gap 
between the ground state and the first excited state. 


IV’fV =afe.n \lk,l-k) + Pk^n |0fc,0_fc) (10) 

with 

gfc.n + 4"^ 

^l,n + 

^ ^ ■ ( 11 ) 
\J^k,n + + E'fc" 

Since the Hamiltonian is the sum of the non-interacting terms 
each one of them is acting on a different Hilbert space, 
and the ground state of the total Hamiltonian is consequently 
a tensor product of all 

= (g) • (12) 
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III. THE SPIN CORRELATIONS EUNCTIONS 

To obtain a generic spin correlation function we can adapt 
the strategy that we used to compute the energy density. Then 
applying Wick’s theorem simplifies the issue further 
since it makes it possible to express any multi-body fermionic 
correlation function in terms of two-body correlation func¬ 
tions. More precisely, it is possible to prove 0 that defining 
for each site j, two fermionic operators, Aj and Bj, via 

Aj = Cj+c^j and = cj - c] , (14) 

any spin correlation function can be written as an ordered 
product of these operators. Hence, due to Wick’s theorem, 
any spin correlation function can be written as a combination 
of one- and two-body expectation values involving only oper¬ 
ators Aj and B^. on the same or different sites. With Eq. (El) 
we obtain 



in) 

The associated energy density ^ is the sum El. ' divided 
by the total number of the spins TV. In the thermodynamic 
limit the energy density becomes 

En, 4 > = -/ y/l — cos((n + 2)k) sin(2(^)(i/c . (13) 

Jo 

According to the general theory of continuous phase transi¬ 
tions at zero temperature ll^ the presence of a quantum criti¬ 
cal point is signaled by the divergence of the second derivative 
of the energy density with respect to the Hamiltonian param¬ 
eter. In Eig. [T]the second derivative of the energy density is 
plotted in dependence of (j) and shows a divergence for the 
value (j) = (j)c = tt/4: independent of n. The singularity is 
ultimately due to the vanishing of the energy gap between the 
ground state and the first excited state at the critical value (Jc 
with the modes k = where j runs from 0 to n -F 1 . 


{AJ) = 0 , 

{BJi = 0 , 

{AiAjJ) = 6i}f. , 

{B,Bk) = -S,k , ( 15 ) 

(BiAk) = Gi^k{n,4>) ■ 

The fact that we have that both {AJ) = {Bj) = 0 and 
{AiAk) = {BiBk) = 0 for i k has several important con¬ 
sequences. In fact, let us consider a spin correlation function 
associated with an operator that is the product of many local 
spin operators, each one acting onto different spins, in which 
trj and/or crj appears an odd number of times on different 
sites. To this operator we may associate a fermionic opera¬ 
tor made by a different number of Aj and Bj operators acting 
onto different spins. Therefore, when we apply the Wick’s 
theorem, we have an expectation value of a single fermionic 
operator and/or an expectation value of two operators of the 
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same kind onto different spins. Hence, taking into account 
Eq. dTSb , such spin correlation functions have to vanish. Con¬ 
sequently, the only correlation function that can be different 
from zero are the ones associated with an operator that is a 
product of local spin operators in which both crj and crj ap¬ 
pear an even number of times. 

To obtain the explicit expression of the non-zero spin cor¬ 
relation functions we need to evaluate ^(n, (p). At first we 
note that, in the thermodynamic limit, the must 

be independent from the choice of i and k but may depend 
on their relative distance r = i — k. With eq. we find 
Gi,k{n, (j)) = Grin, (j)) with 

^ / /N 1 /”^cos(fc(n-|-l-|-r))cos(/)—cosffclr—1)) sin(/i 

Gr(n,d>) = — - , - dk . 

Jo — cos((n -I- 2)fc)sin(2(^) 

( 16 ) 

Solving this integral we find that if r ^ /(n-|-2)-|-l, where / is 
an integer number that runs from —c» to oo, then the Gr (n, 4>) 
vanishes for all values of (jj. This fact, as we show in Sec. [Yl 
plays a fundamental role in the behavior of the entanglement 
property among different spins. 

Obviously, from Eq. (fTSl) and the explicit expressions 
Gi^k{n,4') one can recover all spin correlation functions of 
interest. Here we wish to point out some interesting results 
about some specific ones. 

If one allows for a magnetization along the z direction, i.e. 
(cr|), one finds that it equals Go(n, 0) and, therefore, van¬ 
ishes identically for all possible values of (p and n. Let us 
consider two-body spin correlation functions that can be writ¬ 
ten as (erf cTi^r) n = x,y, z. If y coincide with z the 
correlation function can be written as 

(crfcr|+^) = Go(n, (p) - Grin, p)G-rin, p) . (17) 

Since Grin, p) with r ^ li2 + n) + 1 vanishes we find that 

(«+,) = 0 ( 18 ) 

for all values of n and p. Setting y = x,y the spin correlation 
functions are given by the determinant 


X \ 
i-\-rl 


G-iin,p) G-2in,p) ■■■ G-rin,p) 
G-2in,p) G-iin,p) ■■■ G-r+iin,p) 


G-rin,p) G-r+iin,p) ■ ■ ■ G-iin,p) 

(19) 


Giin,p) G2in,p) ••• Grin,p) 
G2in,p) Giin,p) ■■■Gr-iin,p) 


Grin,p) Gr-iin,p)--- Giin,p) 


( 20 ) 


Numerical evaluations reveal that (erf crf^.^) is non-vanishing 
only when r is an integer multiple of n -f 2, in strong contrast 
to the correlation function (erferlf,^), which is always non¬ 
zero. It changes from negative to positive values in the case 
r varies from odd to even values. This is expected due to the 
anti-ferromagnetic nature of the Ising interaction. 
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Figure 2: (Color online) Behavior of magnetic order parameter in 
the Ising phase p > pc plotted for n = 1, 2,..., 6: red (uppermost 
curve) n — 1', blue n = 2; green n = 3; black n = 4; magenta n = 
5; orange (lowest curve) n = 6 . The dots represent the numerical 
results of my for r going to infinity, Eq. (ED, whereas the curves 
corresponds to the guessed function of the staggered magnetization 
Ea.iT2\. 


IV. THE ORDER PARAMETERS 

As we have seen in Sec. m the behavior of the second 
derivative of the ground state energy density shows a phase 
transition al p = pc = /A for all n. Now we character¬ 
ize the properties of these two phases via the help of the spin 
correlation functions (Sec. lull). 

Let us start from the phase p > pc, i.e. when the system 
is dominated by a two-body anti-ferromagnetic Ising interac¬ 
tion along the y spin direction. Due to the Z 2 symmetry of 
the Hamiltonian ([T]l we cannot compute the staggered magne¬ 
tization by directly applying the definition 
since this gives always a vanishing result. Approaching the 
problem we may first evaluate the value of the magnetiza¬ 
tion with respect to its relation to the long distance correlation 
function along the same spin direction, i.e. 

my = j lim (-l)'-(crfcrf+^) . (21) 

y r —>00 

This can be evaluated via the help of Eq. (l20l i. We have com¬ 
puted for different n numerically the quantity (—l)'^(o'f 
with r up to 200, showing that an increase of the distance r 
results only in a very small variation of my (of a factor less 
than 10“®) for each value of p > j. 

The results that we have obtained for different n and p 
are plotted in Eig. (EJ. This shows the presence of an anti¬ 
ferromagnetic phase along the y-direction for p > pc inde¬ 
pendent of the value of n. However, differently from what 
happens for the second derivative of the density of the ground 
state energy, the staggered magnetization shows a clear de¬ 
pendence on n. Analyzing the numerical data we can con¬ 
clude that the staggered magnetization has the following de¬ 
pendence on p > pc and n: 

mijP^ = (1 —tanip)~'^) ® . 


( 22 ) 
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Figure 3: (Color online) Behavior of the string order parameter for 
(f> < 4>c for n = 1, 3, 5: black (upper curve) n = 1; red (middle 
curve) n = 3; blue (lower curve) n = 5. The dots represent the 
numerical results of the string order parameter Sn given in Eq. ( 124b . 
whereas the curves correspond to the behavior of the string order 
parameter defined in Eq. l l27b . 


From that we can deduce the critical exponent P over n 

Pin) = . (23) 


The fact that the critical exponent /3 depends on n means 
that the class of symmetry to which the models given by the 
Hamiltonian ([T]i belongs depends on n. 

The situation changes drastically when we move in the 
phase below the quantum critical point </) < In this phase 
our Hamiltonian is dominated by the many-body interaction 
terms and ruy drops to zero for any n. It is not straightfor¬ 
ward to find a proper candidate or the role of order parameter 
as it was for the anti-ferromagnetic phase discussed above. 
However, after an heavy numerical analysis we were able to 
obtain a clear picture on the ongoing physics of the system. 

For a system with odd n we can define a string order pa¬ 
rameter as 




= . / lim 


' 2^3 


^n+1 


Oal 


r—n+1 


(24) 


where the operator O = r- 2 {n+i)- Fig-[3]the behav¬ 

ior of this string order parameter Sn for n = 1, 3, 5 is plotted. 
The existence of such a non-vanishing string order parameter 
can be traced back to the presence of diverging localizable en¬ 
tanglement im. This signals the presence of a symmetry 
protected topological order. 

For even n we find that the phase is a nematic one thus we 
can define the following order parameter (since the system is 
translation invariant the quantity is understood to not depend 
on the particular i) 

Bn = {0.,n} = • •<+„) (25) 

As in the staggered magnetic order phase Bn cannot be eval¬ 
uated directly since it vanishes (for any even n the operators 



Figure 4: (Colour online) Behavior of the nematic order parameter 
for (j> < pc for n = 2,4, 6; black (upper curve) n = 2; red (middle 
curve) n = 4; blue (lower curve) n = 6. The dots represent the nu¬ 
merical results of the nematic order parameter B„, Eq. ( 126b . whereas 
the curves correspond to nematic order parameter Eq. ( 127b . 


af or cr^ appear an odd number of times). Again we can cir¬ 
cumvent this problem by defining 

Bn = lim {Oi^nOi+r,n) ■ (26) 

V r—^oo 

In Fig.|4]we plotted the behavior of Bn for n = 2,4, 6. 

Analyzing the numerical data obtained for both defined 
string order parameters, Sn and Bn, we find finally the same 
dependence on n and p, i.e. 

5(") = (l-tan((/.)2)'^ 

= {1 - . (27) 

Summarizing all results we can formulate a general concise 
formula for all order parameters of the whole class of models 
given by the Hamiltonian ([T]i: 

n+2 

Order Parameter = (^1 - " . 

(28) 

Moreover, the existence of a duality, i.e. a transformation that 
brings the order parameters before and after the critical point 
in relation, is thus proven. 

In summary, we find that for both phases we can define 
order parameters that each is ruled by the dominated interac¬ 
tions, i.e. Ising interaction or multi-body cluster interaction. 
In the multi-body cluster phase a strong dependence on the 
size of the cluster n -F 2 is present revealing either a nematic 
phase (even n) or a topologically ordered phase (odd n). 


V. THE ENTANGLEMENT PROPERTIES 

In this section we analyze the entanglement properties be¬ 
tween adjacent spins as well as between a block of spins and 
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the remaining part of the chain. Despite the complexity of the 
class of models under investigation we obtain general results 
showing the relevance of the entanglement features in these 
complex matter systems. 

The object of matter is the reduced density of m spins, 
which is obtained by taking the trace over all remaining spins 
of the ground state. Any such reduced density matrix we can 
decompose by the spin correlation functions 


p(-) = ± 

2m / ^ 
ai, ...,ari 


(^r 


(29) 

where ai runs from 0, x, y, z and denotes the identity. 

The next subsection introduces the concept of different 
types of multipartite entanglement. Then we compute the en¬ 
tanglement properties of adjacent spins and the entanglement 
between a block of spins and the remaining part of the chain. 


In particular, note that the following tripartite mixed state 

P= \4’i)AB{'ilJi\AB <8 |o-i)c(gz|c 

i 

+ \Xi)Ac{Xi\AC^\n)B{n\B 


I 

+ ® |Wi)A(WiU (33) 

i 


with Pi,qi,ri > 0 and = 1 is bi-separable 

though it is not bi-separable with respect to a certain splitting. 
This property and the fact that the convex sum of pure states 
is not unique are the reasons why it is hard to detect genuine 
multipartite entanglement, i.e. a state that cannot be written in 
the above form. Consequently, the entanglement characteri¬ 
zation of multi-partite states needs more than the combination 
of bipartite entanglement criteria . 


A. Definition of hierarchies of multipartite separability 

The quantum separability problem reduces for bipartite en¬ 
tangled systems to the question of whether the state is entan¬ 
gled or not. In the multi-partite case the problem is more in¬ 
volved. First, there exist different hierarchies of separability 
since an n-partite entangled state p may be a convex combi¬ 
nation of pure entangled states with maximally k entangled 
particles. Any tripartite pure state, e.g., can be written as 

|V'fc=3) = ® |<?!>c) 

|V'fc=2) = \4'a) ®\4>Bc)^ \(i>B) ®\(I^Ac) 

or |(/)ab) ® \(t>c) 

| V ^ fc = i ) = \iI))abg ( 30 ) 

where k gives the number of partitions dubbed the k- 

separability. In general a pure state is called fe-separable, 
if and only if it can be written as a tensor product of k factors 
lipi), each of which describes one or several subsystems, i.e. 

IT''"') = IV'i) ® IV' 2 ) d • • • 0 IV”;:) = |V’il/'2 ■ ■ ■'0fc)(31) 

A mixed state p is called fc-separable, if and only if it can be 
decomposed into a mixture of fc-separable pure states 

P = (32) 

i 

where all jT'f) are fc-separable (possibly with respect to differ¬ 
ent fc-partitions) and the pi form a probability distribution. An 
n-partite state (pure or mixed) is called fully separable if and 
only if it is n-separable. It is called genuinely multi-partite 
entangled if and only if it is not bi-separable (2-separable). If 
neither of these is the case, the state is called partially multi¬ 
partite entangled or partially multipartite separable. Note that 
obviously a k = 3-separable state is necessarily also k = 2- 
separable, thus fc-separable states have a nested-convex struc¬ 
ture. 


B. Entanglement properties among adjacent spins 

Let us start by analyzing the case of m adjacent spins 
thus having a maximum distance of r = m — 1. Then all 
spin correlation functions can be expressed by Gr{n, </>) with 
— (m — 1) < r < m — 1. 

Theorem 1. If the number of adjacent spins m is smaller than 
the cluster size, i.e. m < n-\-2, then all k < m-partite entan¬ 
glement vanishes, i.e. the reduced state is k < m-separable. 
If the number of adjacent spins equals the cluster size, i.e. 
k = n -\-2, then there exists a finite range of values of for 
which the reduced density matrix to this set of spins is gen¬ 
uinely n -\- 2-partite entangled (plotted in Fi'g. 13). 

Proof: Let us start with m < n -I- 2. In Sec. uni we have 
computed all spin correlation functions Gr{n,(j)) and found 
that they vanish if ?’ ^ {n 2)1 1. The reduced density ma¬ 

trix pm'^ depends only on a single function, i.e. Gi (n, f). This 
implies that only spin correlation functions that are different 
from zero are the ones along the y-direction. Consequently, 
the reduced matrix pm is a mixture of states being eigenvec¬ 
tors to the single-spin operators ctJ. Applying the following 
local unitary operators 

Uj = exp (-*^crj) : (34) 

(n) 

brings the density matrix pfi of m adjacent spins into a diag¬ 
onal form which obviously is separable. 

In the case where the adjacent spins equals the cluster size, 

(n) 

m = n + 2, the reduced density matrix p „_|_2 depends on 
Gi (n, (j)) and G'_(„+i) (n, f) that corresponds to the spin clus¬ 
ter correlation function (cr^O^ „crfe+„+i)- Again applying the 
above defined local unitary operators Uj to each spin we ob¬ 
tain a reduced density matrix pm'^ that has an X-form 1^ . 
i.e. only entries on both diagonals are nonzero. It has been 
shown that if a density matrix can be written in such an X- 
form, the genuine multipartite entanglement can be exactly 
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Figure 5: (Color online) Dependence of the genuine multipartite con¬ 
currence as function of the weighted interactions (f> for dif¬ 

ferent n that runs from 1 (highest curve) to 12 (lowest curve). Note 
that only for n = 1 genuine tripartite entanglement is non-zero be¬ 
fore and after the critical point and, generally, genuine n + 2-partite 
entanglement decrease with increasing cluster size. 


evaluated by a certain measure, dubbed genuine multipartite 
concurrence introduced in Refs. Thus, by applying 

the above defined local unitaries we find the following expres¬ 
sion for the genuine (n + 2)-partite concurrence for any n and 

max[0, (1 - Gi (n, • (G_(„+i) (n, (/)) + l) 

-1] . (35) 

In Fig. 13 we have plotted the genuine (n + 2)-partite 
concurrence for n = 1,..., 12 and find for certain (j) non-zero 
values. Q.E.D. 



Figure 6: (Color online) Behavior of ^ as function of n where 

the maximization is taken for the genuine multipartite concurrence. 
The red dots are the result of the numerical maximization for any n 
while the black line represent the fit, obtained for large n, in Eq. | I361 . 



n 


Figure 7; (Color online) Dependence of the maximum value of the 
genuine multipartite concurrence on the cluster size 

n + 2. The (red) dots are the results of the numerical maximization 
for any n whereas the (black) line represent the fit obtained for large 
n presented in Eq. ou. 


Looking more carefully at the curves, one observes a 
similar behavior for all n and, except for n = 1, a non-zero 
value of genuine (n -F 2)-partite multipartite is only obtained 
in the Ising phase. Moreover, the genuine (n + 2)-partite 
concurrence is always smaller for bigger cluster sizes. That 
proves that the entanglement in the ground state becomes 
robust against the Ising-interaction (remember the ground 
state of ^ = 0 is a graph state and for ^ ^ a totally 

factorized state 11391 - 1^ ). Consequently, higher cluster sizes 
allow for better properties for running quantum algorithms. 

In Fig. |6] and Fig.|7]a deeper analysis of the entanglement 
properties of the reduced matrix can be found concerning the 
maximal value of the weight (/linax of the two interactions 
which corresponds to the maximal reachable value of genuine 
n + 2-partite concurrence Ggm^^^. Both values show a similar 
dependence that for n > 10 is in good approximation given 


by 


= — , (36) 

“ n 

where the numerical coefficients are obtained by a best fit al¬ 
gorithm. Analogously, the point in which the genuine (n + 2)- 
partite entanglement becomes different from zero depends on 
the inverse of n, plotted in Fig. (0, 



TT 6.2 

2 n 


(37) 


From these equations we immediately reveal an interesting re- 
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Figure 8: (Color online) Dependence of ^ as function of n. 

The (red) dots are the result of the numerical result for specific cluster 
sizes n + 2 whereas the (black) line represents the fit result obtained 
for large n presented in Eq. 03. 


lation between and i.e. 



valid for large n. 

In summary, these cluster-Ising models with different clus¬ 
ter sizes have interesting local entanglement properties. There 
is no bipartite, tripartite,..., n + 1 entanglement, but only for 
large enough values of Ising interaction 0 > one finds local 
entanglement, in particular only genuine n + 2-partite multi¬ 
partite entanglement. In Ref. 114^ the authors computed that 
the maximal value of maximal possible entanglement of two 
adjacent spins in a translation invariant chain was found to 
give a (bipartite) concurrence of (7 = 0.434467. This optimal 
value serves for interpreting entanglement values obtained for 
real physical systems. In the very same manner the maxi¬ 
mal values of the multipartite entanglement quantified by the 
above introduced genuine multipartite entanglement measure 
serves as an reference for real physical system exhibiting clus¬ 
ter and Ising interactions. 


(n) 

Figure 9: (Color online) Here the von Neumann entropy 5m , 
Eq. l l39t . in dependence of a block size m for different cluster sizes 
n is plotted from which we fit the numerical solution of the von Neu¬ 
man entropy, see Eq. 145l l. The value of n runs from 1 lowest (black) 
curve to 8 the highest (pink) curve. 


where iTshannon (a^) is the Shannon entropy 

77shannon(a;) = -xlog 2 (a:) - (1 - a:) log 2 (l - x) , (41) 

and Vj is the imaginary part of the eigenvalues of the matrix 


with 


and 


r = 

nm. 


F' = 

m 


(42) 

/ lip 

n_i ••• 



Hi 

Ho ••• 

n_m+2 





(43) 


nm-2 • • • 

Ho / 



Hr- = 


0 Grin,(j)) 


(44) 


C. Entanglement properties between a block of spins and the 
rest of the chain 

Another important property to analyze in multipartite sys¬ 
tems concerns the entanglement features of a block of m spins 
with the rest of the chain and how it classifies to the holomor- 
phic and anti-holomorphic sectors in conformal field theories. 

For that purpose we have to compute the von Neumann en¬ 
tropy of the reduced density matrix of m spins, 

5i")=Tr(pL"Hog2(p^))). (39) 

Using the methods developed in Ref. ll43ll44l] we find 

= £ i^Sha„non () (40) 

J = 1 ^ 


We have evaluated numerically the von Neumann entropy 
for blocks of length ranging from 2 to 200 spins at the critical 
point (j)c for n that runs from 1 to 10. The obtained values of 
the von Neumann entropy are displayed in Fig.|9] 

Analyzing the numerical data we deduce 

5^^ ~ (0.32 -F 0.18 n) log 2 m + consign) (45) 

The multiplicative constant in front of the logarithmic term is 
known to be related to the central charge of the 1 -F 1 dimen¬ 
sional conformal theory describing the critical behavior of the 
chain via the relation 

C + C 

Sm = log 2 m , (46) 

D 

where c and c are the central charges of the so-called holo- 
morphic and anti-holomorphic sectors of the conformal field 
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theory. Due to the existence of a duality in the system under 
investigation we have that c = c and hence, via Eq. (1451) we 
obtain 

c~ 3-(0.32 + 0.18 n) . (47) 

For two quantum one-dimensional systems to belong to the 
same universality class they need to have the same central 
charge. Since in our case we find a dependence on n, this 
central charge, in addition to the critical exponent /3, Eq. (l23t . 
of the order parameters, Eq. (l28l l. proves that the many-body 
cluster-Ising models fall into different classes with respect to 
their symmetries. 

VI. CONCLUSIONS 

In summary, we have analytically solved, characterized and 
analyzed the properties of a family of models that we named 
n-cluster-Ising models. These are models characterized by 
different cluster sizes (n + 2) and different weighted cluster 
interaction and Ising interaction. We proved that there occurs 
a phase transition exactly when both interactions are equally 
weighted and, interestingly, independent of the cluster size. 

With respect to their symmetries the family of models falls 
into different classes proved via the dependence on n of the 
critical exponent /3 of properly defined order parameters and 
the central charge of the holomorphic and anti-holomorphic 
sectors in conformal field theories. In particular we find that 
the cluster phase has very different orderings for odd or even 
cluster size, namely a topological or a nematic order. Since 
nematic order usually shows up only for non-analytically 
solvable systems these cluster-Ising models may become a 
prototype testing model for exploiting the physical potential 
of nematic ordering of spins. 

In the next step we have investigated how the apparent com¬ 
plexity of the ordering translates to the multipartite entangle¬ 


ment properties shared among spins or block of spins with the 
rest of the system. Surprisingly, exactly all reduced density 
matrices with m adjacent spins smaller than the cluster size 
(= n + 2 adjacent spins) posses no entanglement, whereas the 
reduced density matrices for exactly the cluster size (n + 2) 
possesses genuine n + 2-multipartite entanglement if the Ising 
interaction is strong enough, but not maximal (see Fig. |2). 
This absence of bipartite or n— 1-partite multipartite entangle¬ 
ment is very different from other one dimensional spin mod¬ 
els, i.e. as the Ising one ll^ Id^ll or the XF-model 12^ . 
That computation was possible, because the symmetries of the 
Hamiltonian constrain the state space in the Hilbert-space in 
a non-trivial way enabling even the computation of a measure 
of genuine multipartite entanglement. From the quantum in¬ 
formation perspective these results show that increasing the 
cluster size reduces local entanglement and, herewith, the ro¬ 
bustness of the performance of any quantum algorithm. From 
the perspective of comparison of different condensed matter 
systems the family of models serves as a reference system of 
the possible amount of local genuine multipartite entangle¬ 
ment that can be shared. 

Our family of models can be generalized with respect to 
higher dimensions both in space and degrees of freedom 
(higher spins). These models may become a good testing 
ground for non-trivial spin orderings and serve as a prototype 
for studying the potential of a quantum computer. 
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